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Abstract: Carbon-ion radiotherapy (CIRT) offers unique physical and biological 
advantages over photon radiotherapy. However, some materials and devices in the 
CIRT treatment room become radioactive under bombardment by therapeutic carbon- 
ion beams due to nuclear reactions, thereby leading to possible radiation hazards to 
medical staff and additional and unwanted doses to patients. This study assessed the 
level of induced radioactivity in the treatment room of the Heavy-Ion Medical 
Machine (HIMM) in Wuwei. Monte Carlo simulations using PHITS were performed 
for a conservative case under the conditions of maximum beam energy and intensity 
provided by the HIMM facility. The geometry and configuration of Treatment Room 2 
of the HIMM facility in Wuwei were adopted. We evaluated the activation of air, the 
phantom, and the components of the beamline, such as the primary collimator (PC), 
ridge filter (RF), and multileaf collimator (MLC). For air activation, we calculated the 
medical staff immersion external exposure and inhalation internal exposure caused by 
the corresponding radionuclides. For phantom activation, we estimated the additional 
dose to the patient’s family members owing to secondary photons after treatment. In 
addition, the exemption or non-exemption of the component material activation was 
assessed. 

The results showed that external radiation caused by air activation was the main 
source of the annual effective dose at approximately 0.5 mSv/y. The induced 
radioactivity exposure to family members of a patient after CIRT was approximately 
40 uSv, sufficiently lower than the public dose limit of 1 mSv/a. The induced 
radioactivity of the PC, RF, and MLC was all above the exempt levels after the 
devices were retired, whereas the induced radioactivity of the RS and compensator 
could reach the exempt levels after one patient session. Our study indicated that 
medical staff engaged in CIRT should stay away from the high-dose-rate area of 
induced radioactivity along the beam direction, shorten the residence time in the 
treatment room as much as possible, and store the activated components in isolation 
after the equipment is out of use. Thus, this study provides guidance for accurately 
assessing the level of induced radioactivity in the treatment room for CIRT. 
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1 Introduction 


Compared to photons, heavy charged particles such as carbon ions simultaneously 
provide inverted depth-dose distributions and high relative biological effectiveness 
(RBE) for cancer radiotherapy [1]. These properties make heavy ions more 
destructive to the target volume and less injurious to surrounding healthy tissues. 
Statistics from the Particle Therapy Co-Operative Group (PTCOG) show that by the 
end of 2020, more than 290,000 patients have been treated worldwide with particle 
therapy, 40,000 of whom were treated with carbon ions, accounting for 13.8%, 
indicating that carbon-ion radiotherapy (CIRT) plays an increasingly important role in 
particle therapy. China has also been vigorously developing CIRT, and carbon-ion 
particle therapy systems have begun construction in places such as Fujian and 
Zhejiang [2]. As a result, the issue of radiation protection in CIRT has received more 
attention because high-energy therapeutic carbon-ion beams are necessary for deep- 
seated tumor treatment to increase radioactivity. The induced radioactivity caused by 
the bombardment of therapeutic beams on the materials of the treatment room during 
CIRT gives medical staff and patients unwanted radiation doses [3]. Because induced 
radioactivity is often the main source of radiation dose for radiotherapy-related staff, 
it plays a significant role in radiation protection in the treatment room. 

The physical basis of induced radioactivity is a nuclear reaction, and the induced 
radioactivity can only be produced when the conditions of the nuclear reaction are 
satisfied. For incident neutrons, there are nuclear reactions regardless of energy 
because neutrons are unstable outside the nucleus and are always captured by other 
nuclei to stimulate nuclear reactions. However, for other particles, either their energy 
must reach the reaction threshold for neutron generation or their energy must be 
sufficiently large to cause nuclear fragmentation [4]. The energy range for CIRT is 
generally 80-430 MeV/u. Under such high beam energies, the above conditions for 
the occurrence of a nuclear reaction can be fully satisfied. High-energy therapeutic 
carbon-ion beams produce secondary particles owing to inelastic nuclear interactions 
in the devices of beam delivery systems, patients, etc [5,6]. Most, but not all, of the 
radionuclides produced are short lived. Isotopes that are neutron-deficient are likely to 
decay by positron emission or electron capture (EC), and isotopes with excess 
neutrons are likely to decay by B— emission. Because many of these decays result in 
isotopes in an excited state, gamma rays are subsequently emitted, thereby causing 
exposure to those present in the treatment room [7]. The Shanghai Proton and Heavy 
Ion Center (SPHIC) assessed the occupational exposure of medical staff in clinical 
practice due to the radioactivity induced in patients’ tumors to be approximately 0.508 
mSv [8]. Although the dose level caused by the induced radioactivity is low, this 
unwanted radiation should be minimized owing to stochastic effects. Therefore, it is 
necessary to study the generation mechanism and spatial and temporal distribution 
characteristics of the induced radioactivity in the treatment room comprehensively, 
which is important to understand the radiation hazards correctly, formulate protection 
measures appropriately, and ensure safety during radiotherapy and quality assurance. 
This study conducted a more systematic investigation on the induced radioactivity in 
the treatment room for CIRT using the Monte Carlo (MC) method after consulting 
literature [9]. The Institute of Modern Physics (IMP), Chinese Academy of Sciences, 
developed China’s first medical carbon-ion therapy system, called the Heavy Ion 
Medical Machine (HIMM), in Wuwei. At present, the HIMM facility has entered the 
stage of clinical application at Wuwei Heavy Ion Hospital. In this study, the geometry 
and configuration for the MC simulation were established based on the parameters of 
the horizontal beamline with a passive beam delivery system in Treatment Room 2 of 
the HIMM facility in Wuwei. The specific beam model establishment, parameter 


settings, and calculation formulas are explained in detail in the second section of this 
paper. The third section highlights the simulation and calculation results, from which 
some important conclusions of this study were drawn. This study evaluated the 
induced radioactivity and secondary exposure of medical staff, the public, and 
patients’ families in treatment rooms and provided some radiation protection 
recommendations. In addition, the degree of activation of different materials under 
heavy-ion irradiation was assessed to provide a reference for the disposal of out-of- 
service devices. The aim of this study was to provide guidance for assessing the 
induced radioactivity level in a radiotherapy treatment room for CIRT. 


2 Materials and methods 


2.1 Overview of PHITS 


PHITS is an MC program package used for studying the transport of particles and 
heavy ions. The energies of the heavy ions studied ranged from 10 to 100 GeV/u, and 
the particle types included all the ions (nuclei, nucleons, photons, electrons, muons, 
etc.) [10,11]. Because PHITS can transport ions, this MC program package is 
commonly used in heavy-ion therapy, space radiation, and radiation protection of 
heavy-ion accelerators. Applying the CIRT simulation, PHITS has successfully 
reproduced the dose distribution, including out-of-field doses, and the production of 
secondary neutrons has been verified in comparison with experiments [12-15]. 

The calculation of activation in PHITS must be combined with the DCHAIN-PHITS 
program (referred to as DCHAIN for brevity). DCHAIN is a decay chain analysis 
code for simulating the production, buildup, burnup, and decay of nuclides as a 
function of time in any radiation environment. PHITS provides DCHAIN with a 
neutron flux spectrum and nuclides generated by high-energy nuclear reactions. 
DCHAIN performs decay chain analysis, ultimately outputting time-dependent 
nuclide inventories, decay heat, photon spectra, and effective dose rates [16]. 
Toyohara et al. calculated the radioactivity using the PHITS code and related code 
DCHAIN to analyze the experimental results. The results showed that the ratio of the 
simulated radioactivity to the measured radioactivity was within the range 0.9-1.5 
[17]. Therefore, it is reasonable to use PHITS and DCHAIN to calculate the induced 
radioactivity in this study. 

The [T-Dchain] tally is defined in PHITS to generate the input file of DCHAIN. In 
this simulation, we set the intermittent irradiation mode according to the actual 
situation (8 s for one pulse: 2 s beam-on, 6 s beam-off), and the total irradiation time 
was 20 min. The DCHAIN calculation relies on two primary types of nuclear dates: 
neutron reaction cross-sections and decay dates. The data libraries we chose were 
hybrid data libraries. The neutron reaction cross section library was JENDL/AD- 
2017+FENDL/A-3.0+JENDL-4.0+ENDF/B-VIIl0+JEFF-3.3, and the decay data 
library was JENDL/DDF-2015+ENDF/BVII.0+ ЕМРЕ. In this study, we adopted 
PHITS version 3.28 [18]. 


2.2 MC simulation 


There are four treatment rooms at the HIMM facility in the Wuwei Heavy Ion 
Hospital. The first treatment room contains a horizontal beamline with a spot 
scanning system, and the second treatment room provides both horizontal beam and 
vertical passive beam delivery systems. The third treatment room provides a vertical 
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beamline, and the fourth room contains a 45° oblique beamline, both of which are 
equipped with passive beam delivery systems [19] (Fig. 1(а)). In this study, a 
simplified horizontal beamline from the second treatment room was established as the 
model basis. The dimensions of the treatment room model were 1130 cm х 1698 cm x 
600 cm, and the filling material of the wall was standard concrete with a density of 
2.3 g/em3. The model was determined with reference to computer-aided design 
(CAD) drawings of the treatment room [20]. The simulation of the beamline included 
a pre-collimator (PC), ridge filter (RF), and multileaf collimator (MLC). The phantom 
was a 20 cm х 20 cm х 30 cm water tank irradiated by a 15 cm х 15 cm irradiation 
field with 2-cm SOBP (Spread-out Bragg Peak), which was located at the isocenter of 
the treatment room. Several idiosyncratic components (such as the compensator and 
range shifter) were omitted from the simulations for general calculation. However, all 
of these idiosyncratic components were considered for the evaluation of material 
activation. In conclusion, we adapted two simplified beamlines to achieve the 
different aims of this study (Fig. 1(b) and (c)). The transport of 1 х 1078 carbon ions 
from the surface of the phase space file was simulated to achieve statistical 
uncertainties of within 5% for most statistics, although this came at the cost of 
increasing the CPU computing time [21]. (At some measurement points, fewer 
particles pass through, and even a large number of the simulations cannot reduce the 
measurement error at these points.) For conservative consideration, the maximum 
beam current parameters of the HIMM facility in Wuwei were selected for this 
simulation. A beam energy of 400 MeV/u was chosen, and the current intensity was 4 
х 1058 pps. The main beam parameters of the HIMM are listed in Table 1, where the 
maximum depth corresponds to the range of the carbon-ion beam in water at a 
maximum energy of 400 MeV/u accelerated by the HIMM facility. 
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(b) Simplified horizontal beamline for general calculation 
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(c) Simplified horizontal beamline for material activation assessment 
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(4) Distribution of measurement points 
Fig. 1. MC simulation 


Table 1. Main beam parameters of the HIMM 


Ion Type 
Maximum Energy 400 MeV/u 
Beam Intensity 2 x —4 x pps 
Maximum Depth 27.0 mm 
Dose Uniformity 95% (50 x 50mm’) 
Dose Rate 0.001-1.0 Gy/s 
Beam Size 5:12.0 mm 
Depth Step 2.0 mm 


To achieve the research objectives, some measurement points must be set in the 
simulation. The measurement points were set up every 50 cm from the isocenter in 
three different directions for a total of nine measuring points (the distribution of the 
measurement points is shown in Fig. 1(4)). The three directions were the beam 
direction, the direction deviating from the beam by 45° to the right, and the vertical 
beam direction to the right. In other words, only the radiation field on one side was 
considered. The temporal curves of the dose rates of induced radioactivity at these 
nine measurement points were output to describe the spatial distribution and temporal 
change in the induced radioactivity in the treatment room. In addition, 30 cm from the 
surface of the water phantom was also a point of interest for evaluating the dose 
caused by the activated phantom to the patient’s family members. 


2.3 Data analysis method 


2.3.1 Calculation of the dose caused by activated air to medical staff 


The dose caused by activated air to medical staff includes the internal radiation dose 
(can be divided into ingestion and inhalation) and external radiation dose (can be 
divided into air immersion external irradiation, surface deposition external irradiation, 
and water immersion external irradiation). In this study, combined with the actual 
situation in the treatment room, we considered only inhalation internal irradiation and 
air immersion external irradiation. The activity-dose conversion factor of the air 
immersion external exposure and inhalation dose conversion factor used in the 
calculation process can be found in the ICRP reports [22,23] or GB18871-2002 [24]. 
The dose due to air activation in the treatment room is often characterized by 
estimating the annual effective dose, which is the main source of professional 
radiation doses to medical staff. The external immersion irradiation dose is 
considerably higher than the internal radiation dose. Therefore, the annual effective 
dose to the medical staff from external exposure to air immersion was the focus of 
this study to assess the level of induced radioactivity in the treatment room. The 
radionuclides in the air of the treatment room are mainly derived from thermal 
neutron capture. First, the concentration of each radionuclide induced by radiation in 
air was calculated using the PHITS simulation program. Then, by superimposing the 
exposure dose from each induced radionuclide, the external immersion irradiation 
dose to the medical staff caused by the activated air was obtained. In this study, the 
concentration of each radionuclide in air was compared with the corresponding 
derived air concentration (DAC) limit. 
The calculation formula of the annual effective dose of external exposure caused by 
activated air to the medical staff is as follows: 
, 0) 
where is the annual effective dose of external exposure caused by activated air (in 
units of Sv/a), is the annual working time of the medical staff (in units of h), is the 
activity of radionuclides (in units of Bq/m*), and is the air immersion external 
radiation dose conversion factor (in units of (Sv/s)/(Bq/m°)). Because there was an 
exhaust system in the treatment room, the radionuclide decay constant was corrected 
for in this dynamic situation to 
: (2) 
where is the decay constant of the corresponding radionuclide i (in units of ), is the 
exhaust air velocity of the site, the exhaust rate of Treatment Room 2 is 1.25 /s (4500 
m°*/h), and is the site volume. At the same time, the waiting time (in units of min) for 
the medical staff to enter the room after beam stopping and the stay time t (in units of 
min) for each entry into the treatment room are introduced. Eq. (1) then becomes 
» (3) 
where is the number of times the medical staff enter the treatment room per year, 
which is 5000 per year (20 per day). 
The annual effective dose due to the inhalation of air containing radionuclides (Sv/a) 
can be estimated using the following formula: 


А (4) 
where AAs the inhaled air nuclide concentration (in units of Bq/m*), is the inhaled 


dose conversion factor (in units of Sv/Bq), and is the annual volume of air breathed 


(in units of т/а), calculated at a rate of 20 L/min. This part of the calculation does 


not consider the decay of radionuclides in the air when staying in the treatment room 
and uses the activity concentration of the radionuclides when entering the treatment 
room to obtain a conservative value of the annual effective dose of inhalation internal 
irradiation. 


2.3.2 Calculation of the dose caused by the activated phantom to family members 


To evaluate the dose to the patient’s family members, the following formula was 
applied: 
, (5) 

where is ће gamma-ray dose rate 30 cm from the surface of the water phantom. In 
addition, the following scenario was assumed: the patient leaves the treatment room 2 
min after the end of irradiation, and a member of his/her family attends him/her for 2 
h. Therefore, tı is 2 min, and б is 2 h. Patients receiving CIRT are typically irradiated 
in fractions of 20—30 irradiations, at most. The total exposure was obtained on the safe 
side by multiplying the integrated value of Eq. (5) by 30 [25]. 


2.3.3 Exemption 


If the activity of a radionuclide is less than its exemption value, it is called an exempt 
nuclide. If a substance contains more than one radionuclide and the sum of the ratio of 
the activity of each radionuclide to their respective exemption value is less than one, 
the substance is also exempt. In other words, we must determine whether the 
following equation holds: 
, (6) 

The exempt activity of radionuclides (() can be found in GB1887-2002. In this study, 
the degree of activation of the material was roughly determined by calculating and 
analyzing whether the sum of the ratios satisfied the above equation [26]. 


3 Results and discussion 


In the present study, three aspects of the induced radioactivity in the treatment room 
were evaluated: the activation of air in the treatment room, the phantom, and the 
components of the beamline. After treatment, the medical staff must enter the 
treatment room for further surgery. The activated air in the treatment room causes 
external immersion irradiation to the staff, while the inhaled activated air causes 
internal irradiation. We estimated the occupational exposures of these two 
components for different waiting and lingering times and compared them with the 
national standard. The patient’s family might also be exposed to additional radiation 
due to patient activation after treatment; therefore, we estimated this risk and 
compared it with the standard. The degree of activation of the equipment guides its 
decommissioning; therefore, we also performed an exemption evaluation of the 
beamline components. This study can provide a reference for the treatment and 
disposal of beamline equipment after treatment room dismantling, and we also 
provide our own recommendations. 


3.1 Estimation of air activation 
3.1.1 Overview of the air activation situation 


As mentioned in the second section, we set up nine measurement points to reflect the 
spatial distribution of the dose rate of the induced radiation in the treatment room and 
its variation over time. Figure 2 illustrates the dose rates at each point at different 
times: 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 25, and 30 min after irradiation. As shown 
in these charts, the dose rates of each measurement point in the beam direction were 
significantly higher than the corresponding measurement points in the 45° direction 
and those in the direction perpendicular to the beam direction. The dose rates of the 
latter two orientations were similar, except at the point of 50 cm, where the dose rate 
in the 45° direction was significantly higher than that at other points. The dose rate 
decreased with increasing distance. However, in the 0° direction, the decreasing trend 
of the dose rates with increasing distance was not obvious, unlike in the other two 
directions. These facts coincide with the particle flux distribution in the treatment 
room. 
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Fig. 2. Variation in the dose rate with time at the nine measurement points in the 
treatment room 


The trend of the change in dose rate over time was consistent with the trend of change 
with distance, which decreased with time. Furthermore, the dose rates of each 
measurement point in the 45° and 90° directions decreased sharply within 1 min after 
beam-off and then slowed down. A large number of radionuclides are generated owing 
to nuclear reactions during irradiation, and most of these nuclides are short-lived with 
a half-life of less than 1 min. Owing to the low scattering of heavy ions, there were 
few particles in the non-beam direction, and the probability of a nuclear reaction was 


also low. The effect of short-lived radionuclides in the non-beam direction on the 
induced radiation dose rate was considerably greater than that in the beam direction. 
In the 0° beam direction, the time course of the dose rate was not as fast as in the 
other two directions, as shown in Fig. 2. 

Overall, the distribution of the induced radioactivity in the treatment room was not 
uniform, and the higher dose-rate area of the residual radiation field was mainly 
concentrated along the beam direction. Based on the trends from the simulation data 
and considering the principle of optimization of radiological protection, it is 
recommended that medical staff should try to avoid the high dose-rate area near the 
isocenter along the beam direction when entering the treatment room. 


3.1.2 Annual effective dose to medical staff caused by air activation 


The activity concentrations of each radionuclide at the end of irradiation were first 
compared with the DAC limits. DAC is the limit of the concentration of radionuclides 
in the air based on the annual intake limit (ALI), which indicates that the exposure of 
occupational personnel working in this environment should comply with regulations 
on the dose equivalent limit [26,27]. Compared with the corresponding DAC limits, 
the activity concentration of each nuclide was lower than the corresponding DAC 
values in this simulation. This simulation did not consider the indoor ventilation 
problem during treatment and obtained a relatively conservative nuclide activity 
concentration. In this case, the phenomenon in which the radionuclide concentration 
is lower than the corresponding DAC value can prove that the radionuclide 
distribution in the air of the treatment room after CIRT complies with radiation 
protection regulations. 


Table 2. Radioactivity of main nuclides and their corresponding DAC values 


Nuclide À Activity DAC 
0 Concentration (Sv/s)/(Bq/m’) Sv/Bq (Bq/m?) 
(Bq/m’) 
H-3 2.777810 1.0269x107 / 1.8000х 1077 8.00х107 
Ве-7 2.777910? 3.7382x107 2.7400х 1075 5.2000х 107" 3.00х 107 
С-11 3.3445х10° 2.9122х10° 7.9100х1075 2.2000х1077 1.0010’ 
C-14 2.777810? 1.4547х10* 2.4300x10°'° 6.500010"? / 
N-13 3.9371х10° 2.8093х10° 8.680010" / 7.00х10* 
О-15 8.4500х10° 6.3413x10° 1.0400x10°8 / 7.00х10* 
S-35 2.7779х10° 4.2487х10" 2.9200х10"% 1.1000х10° / 
С1-38 3.0880х10° 2.0354х107 1.9400х1077 7.0300х107" 6.00х10° 
C1-39 2.9856х10° 1.5180x10" 1.3600x1078 7.6000х 107" 8.00х107 
Ar-41 2.8832х10° 7.1297 1.0100х1077 / / 


Table 3. Annual effective dose to medical staff due to radioactivity of шаш nuclides 


activated in indoor ап 


Waiting time Residence time Annual effective dose 

t External radiation Internal radiation 

(min) (min) (mSv/a) (mSv/a) 
0 10 0.3790 1.695x10* 

0 20 0.4077 3.390107 

0 30 0.4107 5.085104 

5 10 0.0948 6.240х10% 

5 20 0.1041 1.248107 

5 30 0.1051 1.872104 
15 10 0.0092 8.500x10° 
15 20 0.0102 1.690х10% 
15 30 0.0104 2.540х10° 


In the second section, we present a formula for calculating the external immersion 
irradiation dose and internal inhalation irradiation dose in the treatment room. 
According to the formula given above, the annual effective doses of external 
immersion irradiation and internal inhalation irradiation caused by air activation in the 
treatment room to the medical staff were calculated, and the results are listed in Table 
3. The number of times the medical staff enter the treatment room was assumed to be 
5000 per year (or 20 per day). The value of the internal inhalation irradiation dose was 
three orders of magnitude smaller than the external radiation dose, which is still under 
the condition of a conservative calculation. Therefore, we believe that the dose 
generated by activated air principally originates from external irradiation. If medical 
staff enter the treatment room immediately after finishing the treatment, the annual 
effective dose of external radiation is 0.379-0.4107 mSv caused by air activation; 5 
and 10 min after beam-off, the annual effective dose is reduced to 0.0948-0.1051 
mSv and 0.0092-0.0104 mSv, respectively. According to Chinese Standard GB18871- 
2002, the occupational annual dose limit is 20 mSv [24]. None of these data exceeded 
the personal dose limits for professionals stipulated in the basic standards. 

Overall, the external exposure dose is an important consideration when calculating the 
annual effective dose owing to air activation. From the estimation results, it was 
concluded that suspending entry into the treatment room after the stop could result in 
a reduction in the annual effective dose of external immersion irradiation caused by 
activated air to medical staff. Therefore, a delay of an appropriate period after 
radiotherapy is recommended by the optimization principle of radiation protection 
before radiotherapy staff enter the treatment room. 


3.2 Estimation of exposure to the patient’s family members from the activated 
phantom 


Patient activation is not traditionally of concern with older photon-based treatments, 
but is becoming increasingly relevant as high-energy ion treatments have gained 


prominence because their spallation products and reactions from secondary neutrons 
result in notable activation [16]. In this study, the patient’s own activation and 
possible additional dose to relatives after receiving CIRT were estimated based on the 
activation of the phantom used in the simulation. 

A large number of radionuclides were generated within the phantom during the 
simulation. Most of them were short-lived nuclides, which existed only during beam 
irradiation before decaying. After therapeutic irradiation, the main nuclides in the 
phantom were O-15, O-14, C-11, N-13, and C-10, which originated mainly from the 
spallation reaction of O-16 caused by high-energy neutrons. After rapid decay of O- 
15, 0-14, and C-10, С-11 and N-13 became the dominant radionuclides 
approximately 10 min after irradiation. In addition, there was also an increase in the 
proportion of other nuclides after 20 min, as shown in Fig. 3. These were mainly 
nuclides with relatively long half-lives, such as H-3 and Be-7 
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Fig. 3. Simulation results of induced radionuclides and the ratio of their activities 
from 0 min to 60 min after irradiation 


Table 4. Main radionuclides in the phantom after irradiation 


Nuclide Half-life (s) Decay Mode 
O-15 1.2210? EC, В" 
O-14 7.06х10! ЕС, В" 
С-11 1.22х10° EC, В" 
N-13 5.98х10° EC, В" 
С-10 1.93х10' ЕС, В" 
Н-3 3.89x108 B 
Be-7 4.60х 109 ЕС 


As shown in Table 4, ЕС and B decay were the dominant decay modes, whereas 
positron annihilation from B+ decay emitted photons. Because beta rays are less 
penetrating, only the secondary photon dose due to material activation was considered 
while considering the water phantom activation. The calculated photon spectra from 
the irradiation phantom were used as the input source term for a secondary-stage 
calculation to assess the dose from the activated materials. The time-dependent 
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changes in the secondary photon dose rate due to phantom material activation in the 
treatment room are shown in Fig. 4. When the phantom was radioactivated by carbon 
ions, it emitted radiation rays into its surroundings from the center of the phantom. 
The secondary photon dose decreased with increasing distance from the phantom. The 
maze structure of the treatment room has a good shielding effect, which significantly 
reduces the dose rates at the entrance and corridor of the maze. However, there is no 
similar shielding structure between family members and patients after CIRT, and the 
photon dose-rate levels around the patient are at a relatively high level immediately 
after one exposure; therefore, it is necessary to estimate the additional dose to family 
members due to activation of patients who receive CIRT. 
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Fig. 4. Variation in the secondary photon dose-rate distribution with time in the 
treatment room 


The exposure to family members was calculated as 40 uSv for the HIMM using Eq. 
(5). Compared with 23.5 uSv for the HIMAC (Heavy Ion Medial Accelerator in 
Chiba) of the National Institute of Radiological Sciences (NIRS) and 20.8 uSv for the 
Hyogo Ion Beam Medical Center (HIBMC) [25], our calculation was almost twice as 
high. This is mostly because we adopted the maximum beam parameters of the 


HIMM in the simulation and obtained relatively conservative data. Most radioactive 
nuclides produced in phantoms have very short half-lives. Therefore, even if family 
members attend to the patient for a prolonged time, the exposure hardly increases. In 
ICRP publication 103, this type of exposure is classified as medical exposure; 
therefore, the public dose limitation is not suitable for use in this case [28]. However, 
based on the results of our calculation, the exposure to the family members of a 
patient is even far less than the public dose limit of 1 mSv/a, therefore this extra 
radiation doesn't have much impact on the patient's family. 
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Fig. 5. Variation in the secondary photon dose rate with time at the measurement point 


3.3 Material activation analysis 


Primary and secondary particles react with the components of the beamline to 
generate large amounts of radionuclides, which activate the material. Assessing the 
activation level of the beamline components is important to guide equipment 
decommissioning. To simplify the simulation, simple geometries were used to replace 
the complex beamline components, and the specific parameters are listed in Table 5. 
We divided the components into two categories: fixed and nonfixed components. For 
the fixed components (the PC, RF, and MLC), we set a continuous exposure of 30 
years using the maximum flow intensity. For range shifters and compensators with 
patient specificity, we only considered the activation after one patient session. 
Continuous irradiation for 20 min was followed by cooling for 24 h. This was 
repeated 30 times for one patient session. 

For the fixed components, a large number of radionuclides were generated after 
irradiation, particularly the MLC. Many of these radionuclides were short-lived. 
However, these are less important for radiation protection. Therefore, considering the 
half-life and activity together, we organized meaningful radionuclides in the table. 
The RS and compensator produced considerably fewer radionuclides. 


Table 5. Material of each equipment in the beamline 


Component Material Density (g/em*)  Ореп field (ст?) 
PC Cu 8.960 15x15 

RF Al 2.700 / 

RS PMMA 1.190 / 


MLC Cu 5% W 95% 18.230 15x15 


Compensator Polythylene 0.942 / 


Table 6. Main nuclides of the PC 
Nuclide Activity (Bq) Exempt Activity/exempt Half-life (s) 


Ni-63 4.06х10° 1.00х10* 4.06x10° 3.19х10° 
Н-3 4.66x10° 1.00х10° 4.66х10° 3.89х10* 
Со-60 9.61х10° 1.00х10° 9.61 1.66х10* 
Fe-55 9.17x10° 1.00x10° 9.17х107 8.66х107 
Ма-22 5.08х10* 1.00х10° 5.08х107 8.21х10” 
Мп-54 8.59x10° 1.00x10° 8.59x10" 2.70х107 
Со-57 1.06х10° 1.00x10° 1.06 2.35х107 
Zn-65 1.29х10° 1.00х 109 1.29х107 2.11х107 
Са-45 8.52х10* 1.00х 10” 8.52x10° 1.41x10’ 
5-35 1.19х107 1.00х10* 1.19х10° 7.56x10° 
Sc-46 2.32x10° 1.00х 109 2.32x10" 7.24x10° 
Co-56 2.61х10° 1.00х 107 2.61 6.67x10° 
Co-58 1.43x10° 1.00x10° 1.43 6.12x10° 
Be-7 1.43х10° 1.00х 10” 1.43x10° 4.60x10° 
Fe-59 1.28x10° 1.00х 109 1.28х107 3.84x10° 
A-137 1.15х10° 1.00х 10" 1.15x10° 3.03x10° 
С-г51 6.08х10° 1.00х10” 6.08х107 2.39x10° 
P-33 1.36х107 1.00х10* 1.36х10° 2.19x10° 
V-48 2.65х10° 1.00x10° 2.65 1.38x10° 
P-32 1.81х10° 1.00х 107 1.81 1.23x10° 
Mn-52 2.12х10° 1.00x10° 2.12 4.83х10° 
Са-47 1.19х10* 1.00х 109 1.19х107 3.92x10° 
Sc-47 1.84х10° 1.00x10° 1.84х107 2.89x10° 
Sc-48 6.80х10* 1.00х 107 6.80х107 1.57х 107 
К-43 6.53х10* 1.00х10° 6.53х10° 8.03 10° 
Co-55 4.80х10* 1.00x10° 4.80х107 6.31х107 
Ма-24 5.71х10* 1.00x10° 5.71x10" 5.40x 10° 
Cu-64 3.60*10° 1.00x10° 3.60 4.57x 10° 
K-42 1.28х10° 1.00х10° 1,.28х107 4.45х10* 
Fe-52 5.19x10° 1.00x10° 5.19х10° 2.98x 10° 
51-31 6.03х10* 1.00х 109 6.03х107 9.44x10° 
Mn-56 2.38х10° 1.00x10° 2.38 9.28x 10? 
Ni-65 1.69х10° 1.00x10° 1.69х107 9.06х107 
Е-18 6.99х10* 1.00x10° 6.99x 107 6.59х10° 
Ат-41 2.35х10* 1.00х 109 2.35x10° 6.58x 10° 
Co-61 5.12x10° 1.00x10° 5,12*10" 5.94х10° 
Mn-51 4.51х10* 1.00х10° 4.51x10" 2.77х10° 
С1-38 5.74х10* 1.00х 107 5.74x107 2.23х10° 
Mn-52m 5.32х10° 1.00x10° 5.32х10° 1.27x10° 
O-15 1.83х10* 1.00х10° 1.83x10° 1.22х10° 


Table 7. Main nuclides of ће MLC 
Nuclide Activity (Bq) | Exempt  Activity/exempt Half-life (s) 


Ni-63 
H-3 
Co-60 
Fe-55 
Na-22 
Mn-54 
Co-57 
Zn-65 
Са-45 
Се-139 
W-181 
Ta-182 
5-35 
5с-46 
Со-56 
W-185 
Co-58 
Be-7 
Fe-59 
Cr-51 
V-18 
P-32 
Mn-52 
5с-47 
5с-48 
W-187 
K-43 
Co-55 
Na-24 
Cu-64 
K-42 
Si-31 
Mn-56 
Ni-65 
F-18 
Co-61 
Mn-51 
С1-38 


2.64х10° 
1.82x10° 
5.81х10° 
4.91х10° 
4.68х10° 
4.48х10° 
6.03х10° 
7.36х10* 
2.91х10* 
1.96x10° 
3.08х10° 
6.47х10* 
2.63х10* 
8.59х10* 
1.37x10° 
2.49x 10° 
8.46х107 
1.80х 10” 
7.43х10* 
2.80х10° 
1.03х10° 
3.19х10* 
9.67х10* 
7.05х10* 
2.63х10* 
2.57x10* 
2.02х10* 
2.27х10* 
5.52х10° 
2.22x10° 
3.79х10* 
1.04х10* 
1.26x10° 
1.0310° 
4.81х10° 
3.29x10° 
1.95х10* 
1.40х10* 


Table 8. Маш nuclides of the RF 


Nuclide 

H-3 
Be-7 
C-14 
O-15 
F-18 

Na-22 

Na-24 


Activity (Bq) 
3.23x10° 
2.99x10° 
2.60* 10° 
2.01х10° 
6.52x10° 
6.92x10° 
7.39x10° 


1.00х10* 
1.00х10° 
1.00х10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00х107 
1.00x10° 
1.00x10’ 
1.00x10* 
1.00х10* 
1.00x10° 
1.00х10° 
1.00х107 
1.00x10° 
1.00х107 
1.00x10° 
1.00х107 
1.00x10° 
1.00х10° 
1.00х10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00x10° 
1.00x10° 


Exempt 
1.00х10° 
1.00х107 
1.00х107 
1.00х10° 
1.00x10° 
1.00x10° 
1.00x10° 


2.64х10° 
1.82x10° 
5.81 
4.91x10" 
4.68х107 
4.48x10" 
6.03х107 
7.36х107 
2.91x10° 
1.96x10° 
3.08х107 
6.47 
2.63х10* 
8.59х107 
1.37 
2.49х107 
8.46х107 
1.80х107 
7.43х107 
2.80х107 
1.03 
3.19х107 
9.67x10" 
7.05х10° 
2.63x10" 
2.5710" 
2.02х10° 
2.27х107 
5.52х10° 
2.22 
3.79х10° 
1.04х107 
1.26 
1.03х107 
4.81x10° 
3.29x10" 
1.9510" 
1.40x10" 


Activity/exempt 
3.23х107 
2.99х102 
2.60х107 
2.01х10* 
6.52x10" 
6.92x10" 

7.39 


3.1910" 
3.89x 108 
1.66x 10° 
8.66% 107 
8.21107 
2.70х107 
2.35х107 
2.11х107 
1.41x10’ 
1.1910’ 
1.05x10’ 
9.91x10° 
7.56x 10° 
7.24x10° 
6.67% 10° 
2.49х10° 
6.12х107 
4.60х 109 
3.84x10° 
2.39x10° 
1.38x10° 
1.23x10° 
4.83х10° 
2.89х10° 
1.57х10° 
8.64 10° 
8.03х10* 
6.31х10/ 
5.40* 10° 
4.57x 10° 
4.45х10* 
9.44x10° 
9.28x10° 
9.06x 10? 
6.59х10° 
5.94х10° 
2.77х10° 
2.23х10° 


Half-life (s) 
3.89108 
4.60х 10° 
1.80х10" 
1.22x10? 
6.59x10° 
8.21x10" 
5.40х 105 


51-31 3.00х10' 1.00x10° 3.00х107 9.44x10° 
P-32 4.99х10! 1.00х10° 4.99x10* 1.23x10° 


Table 9. Main nuclides of the RS 


Nuclide Activity (Bq) Exempt  Activity/exempt ua 
H-3 2.23х10° 1.00х10° 2.23х107 3.89х10* 
Ве-7 4.02x10° 1.00х10” 4.02х10* 4.60х 109 
С-11 1.72х10* / / 1.22x10° 
C-14 1.65x10° 1.0010’ 1.65x10° 1.80х10" 
N-13 4.55х10' / / 5.98х10? 
O-15 231x107 1.00х10° 2.31810" 1.2210? 
F-18 1.3410? 1.00x10° 1.3410 6.59x10° 
Na-22 2.41x10° 1.00x10° 2.41х10° 8.21х10' 
51-31 6.92х1077 1.00х107 6.92х1072 9.44x10° 


Table 10. Main nuclides of the compensator 


Nuclide Activity (Bq) Exempt  Activity/exempt ae 
H-3 7.23х10' 1.00х10° 7.23х10* 3.89х10" 
Ве-7 1.36x10° 1.0010’ 1.36x10* 4.60x10° 
Be-10 3.18х107 / / 4.77x10" 
C-11 6.23х10° / / 1.22х10° 
С-14 6.70х10° 1.00х10” 6.70х 1077 1.80х10" 
N-13 3.3110" / / 5.98х10? 
O-15 7.97х 1075 1.00х10° 7.97x10°° 1.2210? 
F-18 1.11 1.00х10° 1.11x10° 6.59х10° 
Р-33 6.73х10° 1.00х10* 6.73х1077 2.19x10° 


For the PC, MLC, and RF, the sum of the ratios of individual nuclides and their 
corresponding exempt activities was considerably greater than 1. Therefore, the fixed 
components of the beamline after decommissioning the treatment room must be 
stored in isolation and disposed of until their radioactivity levels reach exempt levels. 
The RS and compensator produce weak activity radioactivity and fewer long-lived 
nuclides because they are materials with a small atomic number. Their activation 
levels after a patient session were also exempt. However, considering the As Low As 
Reasonably Achievable (ALARA) principle, we still recommend that they should be 
placed in isolation at the end of treatment. 


4 Conclusion 


The investigation of non-therapeutic radiation is of vital importance for better 
radiotherapy efficacy. In a previous study, we analyzed secondary neutron doses under 
a passive beam delivery system and applied an analytical model to predict the ambient 
dose equivalent [29]. In this study, we performed a relatively comprehensive analysis 
of induced radioactivity, which enriched our understanding of non-therapeutic 
radiation in CIRT. Based on the MC particle transport simulation software PHITS, 
this study evaluated the activation of air, the phantom, and the components of the 
beamline in a CIRT treatment room. For air activation, we calculated the medical staff 
external immersion exposure and internal inhalation exposure caused by the 
corresponding radionuclides. For phantom activation, we estimated the additional 


dose to the patient’s family caused by secondary photons. For component material 

activation, we evaluated the exemption or non-exemption. The following conclusions 

were drawn: 

1. After beam-off, the areas with a high dose rate of induced radioactivity in the 
treatment room were mainly located upstream of the beam direction near the 
phantom. After entering the treatment room, medical staff should avoid staying in 
these areas to reduce radiation hazards. 

2. The annual evidence dose caused by air activation does not exceed the personal 
dose limits for professionals stipulated in the basic standards, and the main source 
is air immersion external radiation. Additionally, after irradiation with 
radiotherapy, an appropriate time delay before entering the treatment room is 
recommended. 

3. The maze structure of the treatment room can have a good shielding effect, 
thereby significantly reducing the dose rate of secondary photons caused by the 
activation of the phantom at the maze entrance and corridor. The dose rate near the 
phantom was high immediately after beam-off and then decreased roughly 
exponentially with time. The possible additional dose to the patient’s family, 
derived from the secondary photons of the present treatment fraction, was 
estimated to be approximately 40 uSv, which doesn't have much impact on the 
patient's family. 

4. Beamline components generate large amounts of radionuclides after nuclear 
reactions with primary and secondary particles. Materials with high atomic 
numbers produce high levels of induced radioactivity and generate long-lived 
nuclides. The induced radioactivity of the PC, RF, and MLC is above the exempt 
levels after the devices are retired, and isolation treatment is required. The RS and 
compensator can reach the exempt levels after one patient session but should still 
be placed separately. 

The additional radiation caused by the induction radioactivity in Treatment Room 2 of 
the HIMM is below legal limits. We made recommendations for radiation protection 
based on the analysis of induction radiation in the treatment room, which are also 
applicable to other therapy centers with CIRT. 
There are several limitations to the simulation setup. A simplified beamline model 
was used, and devices such as the patient couch in the treatment room were neglected. 
Additionally, ventilation in the treatment room during irradiation was not considered. 
Nevertheless, we believe that our study may provide a framework for other induced 
radioactivity studies on the CIRT treatment room. Moreover, in future studies, 
experimental measurements will be performed. 
In summary, our study provides references for assessing the induced radioactivity 
levels in treatment rooms for CIRT and provides several radiation protection 
recommendations for medical staff and equipment decommissioning. With the 
popularization of CIRT in China, this study might help different CIRT centers to take 
action when considering the radiation protection of induced radioactivity. For some 
new radiotherapy techniques, such as flash, multi-ion combined therapy, and arc 
therapy, it is also important to re-evaluate the induced radioactivity level under these 
new radiotherapy techniques. The research methods in this paper also provide 
guidance for the study of induced radioactivity under these new techniques so that 
they can be better applied in clinical practice. 
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